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ABSTRACT 


He discuss in detail a model which makes definite predictions 

for the fractionation of isotopes in sputtered material. The 

fractionation patterns can be non-linear , and the pattern for a 

particular set of isotopes depends on the chemical matrix within 

which those isotopes are contained. Calculations are presented 

for all non -mono iso topic elements contained in the minerals perov- 

akite, anorthite, ackermanite, enstatite, and troilite. All isotopes 

are fractionated at the level of approximately 4-6 °/oo per atomic 

mass unit. 0 is always positively fractionated (heavier isotopes 

sputtered preferentially), and heavier elements are generally 

negatively fractionated (lighter isotopes sputtered preferentially). 

The value of 6(*®0:^0) is always less by about 1.8 °/oo than a 

17 16 

linear extrapolation based upon the calculated 6( 0: 0) value 

would suggest. The phenomenon of both negative and positive frac- 
tionation patterns from a single target mineral can be used to make 
an experimental test of the proposed model. 
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Recent experiments [Russell et >1 . , 1980] and theoretical models 
[ Watson and Haff , 1980; Watson, 1980] have addressed the question of 
Isotopic fractionation of surfaces which have been exposed to fluxes of 
low-energy (keV) ions. Earlier work has tended to center mainly on chemical 
composition changes produced by ion bombardment [ Liau et al . . 1977] since 
several convenient techniques (e.g. , Rutherford backs cattering) exist to 
analyse the near-surface composition as a function of depth. Theoretical 
studies of these processes, e.g., by Haff and Switkowski [1976], are always 
faced with the necessity of having to prescribe quantitatively the differ- 
ences in chemical binding energies between atomic species in order to pre- 
dict the degree of chemical fractionation. In experiments involving ion 
sputtering of isotopic components of a single element, the experimenter is 
faced with the complexities imposed by the small size of the expected 
effects, but the theoretical analysis of the sputtering process and the 
determination of relative sputtering yields for a given suite of isotopes 
is made correspondingly easier. 

Based upon an original motivation to construct a sputtering theory of 
binary and more complicated materials, which could be tested against experi- 
ment, Watson [I960] and Watson and Haff [i960] developed a model which de- 
scribed the sharing of recoil energy amongst the components of a given target 
The results of the investigation indicated that little fractionation was 
likely to occur if fractionation effects were due entirely to differences 
in the bulk recoil fluxes of constituent target atoms. The calculated 
fractionation, 8^ ^ , due to non-stoichome tries in the internal fluxes 
alone was only ~ -1 °/oo (parts per thousand) for ^Ca with respect to ^Ca 
in a mineral target. We here define 
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where is the number of target atoms of type i sputtered per incident 
projectile, and is the abundance of atoms of type i at the surface of 
the unsputtered target. Thus we compare the composition of the sputtered 
material to the undisturbed target composition. However, by making the 
critical assumption that the atoms found in the extreme outer layer cannot 
participate in the collision cascade on an equal footing with the internally 
recoiling atoms, Watson [1980] arrived at an expression for the isotopic 
fractionation expected in the material sputtered by the collisional-recoil 
process from a surface containing isotopes of species k and I, 


where 
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n^ is the fractional abundance of isotope i, is the low-energy 
collision cross section between atoms of type i and type k, and 7^ is a 
function of the atomic masses: 
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The sums are taken over all target atomic species, but k and i refer only 
to isotopes of a single chemical element. 

We emphasize that (2) applies only to the material actually sputtered 
away from the target, and not to the composition of the target surface 


This expression for 6. ( is the one given in Watson [1980]. It differs 

K,l 

slightly from that found in Watson and Haff [1980], although numerical 
values are similar in the two cases. For equal cross sections, the two ex- 
pressions are identical. 
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subsequent to sputtering. Modification of the surface composition can be 
described by models which Incorporate both preferential sputtering effects 
and subsurface diffusion processes, but such a project is beyond the scope 
of the work discussed here. The abundance factors n^ appearing in (2) 
refer to the Instantaneous atomic abundances. In general these are not 
constant in time if 6 f 0, and thus 6 ■ B(t). The results which are re- 
ported here , therefore, refer only to low dose experiments where a limited 
amount of material is sputtered from the sample. 

To illustrate the results one would expect if fractionation occurs 
according to (2), we first specialise to several Idealized cases. Consider 
a target composed of a single element, which iu turn is composed of only 
two isotopic species, 1 and 2. The total cross sections are all equal 
to a common value, hence 



n l * n 2 y 12 
°1 7 12 + n 2 


1 . 


(5) 


We let the masses be and Mg ■ + &M for species 1 and 2 respectively, 

and furthermore define AM/M^ * €. To illustrate the fractionation behavior 
for small mass differences we take € « 1. Then g » 1 • J £ 

*1,2 "i* 2 (n l " n 2 } * (6) 


For such a two component system, the sign of the fractionation depends upon 
the abundance factors n^, and not upon the mass values. Thus if n^ > ng in 
the surface layer, 6^g > 0 and species 1 is sputtered preferentially. More- 
over, the fractionation effect is quadratic in the mass difference e. For 
a mass increment of c * 0.1, (6) gives a maximum limiting value of 2.5 °/oo 
for the magnitude of the fractionation effect (see Fig. l). 
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The expression given by (2), and hence the result (6), arise from a 


detailed solution of the transport equation describing energy sharing 
amongst recoiling atoms [ Watson , I960, Watson and Haff . 1980]. A source 
of particles Is created by primary collisions between the incident Ion and 
a target atom. These recoiling atoms are typically of very low energy, 
and interact much in the manner of hard spheres. For projectiles in the 
energy range of a few keV to a few hundred keV, a description of the 
sputtering process based upon such a picture [ Sigmund , 19P9] is well- 
established. Absolute sputtering yields (number of atoms ejected per 
incident ion) of many materials can be computed to an accuracy within a 
factor of two or three, and, more importantly, the dependence of the sput- 
tering yield on incident ion energy, mass, and charge can be reproduced 

adequately by the cascade theory, for a wide variation in the parameters. 

_2 

Furthermore, the predicted E dependence of the yield on the energy of 
the sputtered particle has been independently verified by several investi- 
gators [ Thompson , 1968; Weller and Tombrello , 1978]. 

These calculations have in each case been performed for a uniform 
distribution of scattering centers. In an actual sputtering experiment, 
however, the distribution of atoms fills only a half-space. The actual 
calculations are thus perturbation theory calculations, with the target 
surface introduced only at the final step as the boundary across which the 
sputtering flux should be taken. Sequential collision events such as shown 
in Fig. 2 are Included in the calculation, but cannot physically occur in 
the true target since the collision point is outside the target surface. 
Watson [1980] and Watson and Haff [1980] made an attempt to Include effects 
Introduced by the presence of a surface in a less rigorous way. They 
postulated that the atoms comprising the extreme outer surface layer of the 
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target cannot participate fully in the recoil cascade, primarily because 
the geometry makes it difficult for them to transfer energy to subsequent 
atoms after they themselves have been struck. Sputtered particles are then 
imagined to derive from this passive surface layer as it absorbs energy 
from recoiling atoms deeper in the target. It is this model, which attempts 
to take account of the effects of a non-uniform (i.e., half-space) distribu- 
tion of target atoms, that leads to the expression (2) for the fractionation. 

The "surface-flux" model was adopted because the non-stolchiometrlc 
emission of isotope species expected from the bulk recoil flux alone 
was found to be much smaller than indicated by experiment. Thus Russell 
et al . [1980] found 6(^Ca:^Ca) for material sputtered from a plagioclase 
target to be on the order of 20°/ oo, while predictions based on non- 

stoichiometries in the bulk recoil flux indicated values no less than 
-l^oo [ Watson and Haff , I960]. On the other hand, the surface-flux model 
is in adequate agreement with the data of Russell et al . [1980] for Ca frac- 
tionation in plagioclase and in fluorite. However, the peculiar fractiona- 
tion patterns inherent in the model, as shown clearly by (6) which predicts 
(a) that the value and even the sign of 6 depend upon the abundance factors 
of the target components and (b) that it is possible for the fractionation 
pattern to contain no linear term in the mass Increment, have not yet been 
adequately tested experimentally. In the subsequent discussion we will 
clarify the physics behind the fractionation mechanism at issue here and 
suggest some further experiments designed explicitly to test the model. 

We continue for a moment our discussion of a two component medium. Accord- 
ing to Watson and Haff [1930], the bulk flux of recoiling atoms is very nearly 
stoichiometric (except for extreme differences in target atom masses, or 
for the very highest energy particles which, however, constitute only a small 


portion of the total flux). This flux has an enargy dependence of the fora 

o 

An^/E , for species 1, where A Is a constant. For simplicity, Imagine head- 
on collisions of these bulk-flux components with stationary surface atoms. 
The surface flux of type-1 atoms Is then 


n l °2 
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where the first term describes atoms of type-1 ejected by collision with 

atoms of type-1, and the second term describes atoms of type-1 ejected by 
collision with atoms of type-2. Equation (?) leads to the partial sputtering 
yields 

dSj « nj (nj + 7 12 n 2 ) dE/E 2 (8) 


and similarly 

2 

dS 2 “ n 2 ^ n 2 + 7 21 ”l^ dE / E ' 


(9) 


The fractionation of the sputtered material with respect to the bulk is 
therefore 


. K/ftS; 

1,2 “ n~/n 2 


( 10 ) 


which reduces to (5) upon substition from (8) and (9). The 7 -factors in ( 7 ) 
come from the kinematic limits on the maximum amount of energy which can be 
transferred in an elastic collision. The source of the fractionation in 
this model is therefore seen to arise from the energy transfer mismatch 
between the surface species and the bulk recoil flux; l.e.. It Is easier 
to transfer energy to a similar mass in a collision than to a much different 
mass. The effect is larger in the surface layer than Internally because of 
the fact that surface atoms are allowed to Interact with the cascade in one 
step only: there Is no opportunity for the effect to be averaged away over 
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many collision*. Equations similar to (10) and (5), and their generalisa- 
tion to polyatomic media, were derived by Watson [1980] and Watson and 
Haff [i960] without the simplifying assumption of head-on collisions used 
here. In the general case in which different chemical species are present 
in the target, the relative total collision cross sections which describe 
scattering of the various distinct pairs of atoms also appear in the ex- 
pression for &, as in (2). 

These remarks apply to isotopic components of a single element. If 
more than one chemical species la present, fractionation of one chemical 
species from another will in general occur, but its magnitude is determined 
more by details of the target chemistry than the atomic masses and abundances . 
However, expression (2) for the fractionation applies to any target, regard- 
less of its chemical composition, as long as it is only the fractionation 
amongst isotopic components of a single chemical element that is desired. 
Moreover, since members of any isotopic suite are ejected from the surface 
by means of collisions with all the atomic species in the target, it is clear 
that the fractionation pattern for these isotopes can depend strongly on the 
mass of each kind of target atom. It is this feature which we feel will 
provide the most convincing test of the validity of the fractionation model 
described here. 

As an illustration of the effect of this "background" mass upon the 
fractionation pattern of a particular pair of Isotopes 1 and 2, with masses 
Mj and (Mj < M^, say), we consider a target composed principally of type-3 
background atoms, with abundance n^, but which contains a small proportion 
of the Isotopic species 1 and 2 so that nj+n 2 «nj. Type-3 atoms need not 
be Isotopes of the same element as type-1 and type-2 atoms. The last 
inequality is not necessary for analysis but it leads to a clear picture 
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of the fractionation procaaa. Esaentlaliy all atoms which are ejected will 
arise froa collisions with type -3 atoms. Thus when Mg < Mj, we expect 
species 1 to be preferentially sputtered, so that B g l < 0, but If Mg < Hj, 
then Bg t > 0. For the case M i < < Mg, the fractionation will be small, 

and it will vanish for soma value of Mj in this range. For this case we 
find from (2) 


since the cross sections and •** equal. The fractionation in this 
case Is Independent of the abundances of the Isotopes In question, as long 
as they are much less than unity. If Mg ■ Mj + AM, and c » AM/m^, then for 


c « 1, 


( 1 - "A N 

!<1 ' ll.M,/H,/ ' 


so that the fractionation is linear, with - c <6 . . s e. Figure 3 shows 
Bg j as a function of Mg for the particular choice € » 0 . 1, Mj - 40 (and 
Mg - kk). Also shown is Bg ^ evaluated according to (11), so that an expan- 
sion in powers of c is avoided. The two curves differ significantly in 
the region Mj ■» Mj, but they both illustrate the fact that the magnitude and 

sign of 6„ . can be drastically affected by the choice of the chemical matrix 

*•> 1 

in which the isotopic species of interest are imbedded. If for simplicity 
we take n^, ng, and n^ to refer to Isotopes of a single element, then in 
the special (but unlikely) case Oj » n^ * n^ - l/3 we find Bg i" ^ * n ^ 

Bj^ 1*0 for Mg « Mj + 1 and Mj * Mj + 2 with € * l/Mj « 1. Here the middle- 
mass isotope Mg is better coupled to the underlying recoil cascade than 
either of the end members. 
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Although, with the exception of epee lei exact (tee (6) end the ebove 
paragraph) , 6 la proportional to the mete lncrenent « for small e, If this 
quantity la not small, then 6 will depend upon the increment In e more coopli- 
ceted wey. We do not need to eppeel to en expene ion in c to eee this effect, 
since, within the model (2) Is exect. Figure k gives en illustretlon of non- 
linear it lea In the fractionation, due to violation of the condition c « 1, 

hO k2 

for a pure Ca target containing the terreacrlal abundances of Ca, Ca, 

^Ca, ^Ca, 4 ®Ca, end **®Ca. The curve was calculated from (2). The 6- 

hO 

values with respect to Ca are all negative since the dominant abundance 
hO 

of Ca leads to preferential ejection of the lighter isotopes. 

Ve turn now to more interesting targets, taking as a first example SiOg. 
Because for the light elements stable light isotopes are generally more 
abundant than heavier Isotopes, ve would expect for a pure Si target (n^ ■ 
0.9217, n g9 • 0.0h71, n^ • 0.0312) that 6 2g ^ and Bjq gQ would be negative, 
while for a pure (solid) 0 target (n^g ■ 0.99759, * 0.00037, n^ g • 

0.0020h) and 6^ would also be negative. However, for the SiOg 

target, the qualitative conclusions reached for the sign of the 6-values 
above are no longer all necessarily true. First, although the addition 
of 0 to a Si target would tend to make the 6g g ^ and 3^ gg values even 
store negative, the addition of Si to an 0 target would have the opposite 
effect on jg and 6^ jg, and could even make them positive. There is 
the additional complication of the cross-section values to be used [see (3)}. 
When two or more chemically distinct species are present in the target with 
non-negliglble abundances, the total low-energy scattering cross sections 
help to determine the fractionation values. The enter the calculation 
because the collision probabilities determining the coupling of the surface 
layer to the bulk depend on the product of abundance times cross section. 
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Because of the structure of (5), one of these cross sections any be chosen 
arbitrarily, i.e., we are Interested in the relative, not absolute, sputter- 
ing rates. Consequently, the frsetionetion values depend only upon the 
ratios of cross sections, and not on their nstnitudes. When a target is 
composed of isotopes of one elenent only, or where collisions between such 
isotopes and one other distinct elenent are the only collisions which are 
ioportsnt, then the ratios of all relevant cross sections becone unity, as 
is the case in the above exenples, and & then depends on the masses and 
abundances alone. 


To begin with, we Ulustrste the frsetionetion expected for both SI and 


0 in S10 £ , with the values (taken froc n Born-Mayer nodal of the atom) 

"si-61 <■ "28,28 ■ "28,29' ,tc) ‘ S ‘ 91 **' "si-0 "28,16 * "28,17 

°29 16' * 5. 14 A 9 , and Op ^ p » 6.66 A . V. find froa (2) 



which is essentially a linear relationship. The important point though Is 
that the trend of the fractions*.? on corresponds to preferential emission 
of heavy isotopes for one chemical element (0), and to preferential emission 
of light Isotopes for another (Si). 
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Th« above value* ara dependent upon tha cross aaedona for scattering 
of Si and l» a coas. To gat an idaa of tha sensitivity of & to tha cholca 
of croaa taction, ve arbitrarily intarchanga tha Si ■* Si teat taring croat 
taction with tha 0 -* 0 teat taring croaa taction in tha evaluation of 6 , 
!•••> X , * 5. 81 A j "it 4 0 ■ 5. lfc A at before. 


Then 



*17,16 ' 

3.9 °/oo 

and 




*18,16 ‘ 

5. 9 °/oo, 

with roughly the same 

amount of non-linearity as above, while 


Q ft 

29,28 

-6.9 °/oo 

and 




s 

o 

-ll*.l °/oo, 

still almost linear. 

Finally, in the 

case that all cross sections 

equal to eech othe<:. 

we find 



*17,16 

1* . 1 °/oo , 


*18,16 

6.1* °/oo, 


*29,28 

-6.5 °/oo. 

and 




So. 28 

-13.2 °/oo. 


Figure 5 summarises the fractionation patterns for theae three choices of 
the cross sections. 

The dependence of & upon the (generally not well 'known) collision 
cross sections Is a fact one must live with. It makes a reliable calculation 
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of expected 5-values more difficult. Still, the patterns illustrated in 
Fig. 5 are relatively stable against variations in the o^. One reason is 
that only the relative values of the cross section are required in (3). 

In order to illustrate the application of the above results to partic- 
ular targets, we give some examples of predictions based upon (2) as applied 
to selected high temperature condensate minerals. Mineral targets were 
chosen for several reasons. The most complete characterization of isotopic 
fractionation effects induced by sputtering was recently carried out on 
plagioclase and fluorite targets [ Russell et al . , 1980] . The present model 
was developed in part to try to understand these results. In earlier work 
the question of isotopic fractionation of mineral components during sputter- 
ing was addressed by Switkowski et al . [1977] to try to account for the 
observed abundance of Si and 0 isotopes 5n lunar fines. More recently 
Clayton [1980] has raised the question of whether sputtering of interstellar 
grains could be the source of certain isotopic anomalies observed in meteor- 
itic inclusions. In the examples illustrated below we do not endeavor to 
come to grips with the problem of the likeliness of sputtering as a signifi- 
cant isotopic fractionating agent in pre-solar system history, but it is 
nevertheless interesting to see how appropriate mineral phases [ Grossman , 
1972] would be expected to fractionate in the present model. By presenting 
these examples, we also gain a clearer idea of the magnitude of the effects 
which might be expected in an actual experimental test of the model. 

Figures 6-10 show the predicted fractionation patterns for all sets of 
isotopes in, respectively, perovskite (CaTiO^), plagioclase (anorthite) 
(CaAlgSigOg) f melilite (akermanite (CagMgSi^)), enstatite (MgSi0 3 ), and 
troilite (FeS). To be definite we have set all cross sections 
equal. The patterns exhibit several interesting features. Perhaps 


most striking is the fact that the fractionation effect in one mineral looks 

essential the same as in any other mineral containing the same kinds of 

atoms. The same basic pattern if also encountered in a number of other 

minerals not illustrated here: corundum (AlgO^) , magnetite (Fe^O^), rutile 

(TiOg), eskolaite (Cr g Oj) , various spinels (FeAlgC^, ZnAlgO^, HnAl^Oj. , MgAlgO^), 

gehlenite (CagAlgSiO.^), hibonlte (CaAl^O^g) , fluorite (CaFg) , albite 

(NaAlSi^Op), diops id e (CaMgSigOg), and forsterite (MggSiO^). In each case 

0 is fractionated to a level of roughly 4-6 °/oo per mass unit, with a slight 

(~ 1.6 °/oo) , but remarkably constant, non-linearity, to be discussed below. 

The 0 curves correspond to positive fractionation, and the reason is the same 

16 

as for the 0 component of quarts, Fig. 15. In each mineral 0 is the light- 
est component, and the coupling to the high mass components leads to prefer- 
ential emission of the heavy isotopes. The effect is not a strong function 
of the mass of the heavy atoms, as long as they are reasonably abundant. In 
the O-containing minerals illustrated above, 0 accounts for, to within a few 
percentage points, 60$ of all atoms by number, whilo the heavier atoms, 
usually some combination of Al, Mg, Si, and Ca, make up the remaining 40$. 

Thus 0 fractionation does not vary much from one mineral to the next. For 
similar reasons, the heavier conponents of the mineral are in each case 
negatively fractionated, i.e., the lighter isotopes are sputtered preferen- 
tially. 

The same kind of results are predicted for the non-oxygen containing 
mineral troillte (FeS). S is positively fractionated and Fe negatively 
fractionated, in both cases at the level of a few parts per thousand (ppt) 
per unit mass. 

These seem to be the two hallmarks of sputter-induced isotopic frac- 
tionation in many minerals: (l) all isotopes arc fractionated at the level 
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of a few ppt per unit mass, independently of the precise mineralogy of the 
target> and (2) sputtered 0 is isotopically heavy with respect to the 0 in 
the target! while heavier sputtered elements (Mg, Si, Ca, Ti) are isotopically 
light. The magnitudes of the fractionations predicted are large compared to 
the precision of the most careful mass spectrooetric analyses (—0.1 °/oo 
[Lee, 1979)). The fact that sputtering of Isotopes of the light and heavy 
elemental components of a given target is predicted to produce fractionations 
of different sign sharply distinguishes the kind of process described here 
from mechanisms based on thermal or diffusive kinetics, which can also lead to 
isotopic fractionation. Fractionation due to mechanisms such as these must 
be either positive or negative, but not both, for a given set of processed 
isotopes. The "mixed" positive and negative fractionation patterns produced 
in this model would be suggestive signatures for the origin of any material 
which exhibited them. 

Some of the illustrated fractionation curves exhibit, in addition, 

definite non-linearities. The most regular and for the most part the 

18 

largest of these is shown by the 0 series. In every case studied, 0 is 

less abundant in the sputtered material than a linear relation based on the 

6(17,16) value would suggest. Alternatively, based on the 6(18,16) value, 

17 

we could say that the 0 abundance is enhanced. Our procedure has been to 
compare all 6-values predicted according to (2) with a straight-line value 
obtained by passing a line through the origin and through the B-value for 
the lightest pair (M + A,M) of stable isotopes. The upper panels in Figs. 

5-10 show the non-linearities calculated by this scheme, 

l M+ i e i 6(M+ A,M) - 6(M+ i,M) . (15) 

The ^-values for 6(18,16) are confined to the narrow range +1.8 to + 1.9 °/oo. 
and are essentially independent of the mineralogical matrix. Small non- 
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linear features also appear in some of the heavier isotopic series, but they 
are neither so regular nor generally so large as those found for 0. 

Before pursuing the possibilities inherent in a mechanism which can 
lead to anomalous isotopic patterns, however, we need to establish the valid- 
ity of the model proposed here, and this can be done only by appeal to appro- 
priate experiments. The present model arose out of the Ca fractionation 
data of Russell et al . [i960], Watson and Haff [1980] found it necessary 
to invoke the surface flux model in order to obtain fractionation values of 
the magnitude required. But the model did not predict the experimental re- 
sults; it can be said only to be consistent with them. However, the predicted 
positive-negative fractionation pattern illustrated in the above figure pro- 
vides the kind of yes-no test that can give substantial support for or 
evidence against the proposed sputtering mechanism. 

A test that seems to have merit involves measuring isotope ratios from 
the same chemical element in two different targets. The "background" masses 
of the auxiliary partners would be different in the two cases, and chosen to 
yield an expected positive fractionation of the isotopes of the given element 
in one case and a negative fractionation in the other. Because experience 
has already been gained in the measurement of fractionation of Ca-contalning 
minerals under sputtering conditions [ Russell et al . , 1980] , we examine 
predictions of (2) for two possible Ca-containing targets. 

Figure 11 shows expected fractionations for fluorite, CaFg, and for Cal^. 
The elements F and I are each composed of a single stable isotope with mass 19 
and 127, respectively, and these two masses bracket the range of the Ca iso- 
topes. The Ca sputtered from CaFg is predicted to be strongly fractionated 
in a negative sense, and the Ca sputtered from Cal^ in a positive sense. 

The magnitudes of the 5-values are large, ranging from — 10 °/oo to — 50 °/oo. 
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Moreover, especially in the case o£ Cal g , the non-linearities are also large, 
the maximum 5 -value exceeding 7 °/oo. The direction of the non-linearity is 
the same (positive) for both CaF g and Cal g . This is also the case for the 
Ca-metal fractionation curve illustrated in Fig. 4. A positive §-value for 
a given isotope means that the isotope is more depleted in the sputtered 
material than a straight-line law based on the two lightest isotopes would 
suggest. Thus for Cal,-,, the enrichment of heavy isotopes is less than the 
value of 5(42,40) suggests, and for CaF g , the depletion of heavy isotopes 
exceeds a prediction based upon 6(42,40). This target pair could provide, 
then, an economical test of the sputtering model proposed here, since both 
the positive-negative fractionation feature as well as the non-linearity 
characteristic can be investigated at the same time. 

We conclude the presentation of this sputter-induced fractionation 
process by discussing some of the qualifications and difficulties which 
attend it. First, it must be emphasized again that the 6-fractionation 
values apply strictly to the material which is sputtered away in the initial 
stages of bombardment, not to the surface material which is left behind. 
Furthermore, since 6 depends upon the abundance of the various constituents 
at the surface [the n^ in (3)], the amount of fractionation will change as 
sputtering proceeds. In principle, the n^' s should be considered to be 
functions of time, n^ • n^(t), and therefore 6* 6(t). The total effective 
6 would then be obtained by integration. Unfortunately this is a complicated 
and uncertain procedure. The time dependence of n^(t) depends not only on 
the instantaneous partial sputtering yields of the various species, but it 
depends also upon how the material at the extreme outer surface is mixed 
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with the underlying materiel. This mixing process t haff and Switkowskl . 197?] 


alwaya accompanies Ion bombardment in the region of a concentration gradient, 
and at the present time we are unable to treat the effect in a suitably 
quantitative manner. A rule of thumb which has been found useful is that 
preferential sputtering effects persist, at decreasing levels, until a thick- 
ness of material on the order of the range of the Incident ion has been 
removed. For typical sputtering experiments this distance is quite small, 
on the order of 1700 A in the experiments of Russell et al . [1980]. These 
authors found in fact that 6-values approached aero once 15-50# of the ion 
range had been sputtered away. The 5-values quoted in this paper refer to 

the material removed in the very first moments of sputtering. Subsequent 
sputtering can be expected to leal to a dilution of the effect, and hence 
to smaller effective 6-values for the sputtered material. (Clearly &-*■(> 
rigorously when the target has been entirely sputtered away.) For this 
reason the 6-values predicted here will tend to overestimate the correspond- 
ing measured quantity. Watson and Haff [1980] found that their predicted 
fractionation effects exceeded the measured values typically by a factor on 
the order of 2. 

It is also important to keep in mind that the surface flux model, on 
which all the above results are based, is an idealized and perhaps not 
totally consistent treatment of the effect which the introduction of a half- 
space-type boundary can be expected to have an emission of particles from 
the surface. Thus, if the true composition of the sputtered flux contains, 
in addition to the surface flux, a component arising directly from the 
(essentially stoichiometric) internal recoil cascade, then the predicted 
fractionation effects will be reduced. 

Nevertheless, one can pose a rather definite test of these ideas, as 
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exemplified by the CaF^, Calg system. A negative result in such an experi- 
ment would force a close re-examination of the surface flux model. A 
positive result would provide a significant stimulus to further investiga- 
tions into the fractionation process. Especially interesting are problems 
which need to be addressed about the role of sputtering processes in space. 
Observed isotopic patterns in meteoritic inclusions provide Important clues 
to the origin of the sun and planets. Since sputtering of grains by shock 
waves in the Interstellar medium is likely to have occurred [ Dwek and Scalo , 
1979], we may ask how the resulting fractionations combine with those in- 
duced by nuclear reactions. The answer to such questions will depend upon 
the outcome of experiments designed to test fractionation theories such as 
presented here. 
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FIGURE CAPTIONS 


Fig. 1 


Fig. 2 


Fig. 3 


The 6-value of (6) la plotted for a binary Medium of two iaotopic 
apeciea M^, Mg with abundances n^, r.^. The masses are related 
according to Mg ■ Mj ♦ AM, and c m AJ0lj is taken to have a value 
of 0.1. 

. This schematic picture shows the type of unphysical collision 
(taking place outside the solid) which is unavoidably included in 
most calculations of the sputtering yield. The surface flux model 
attempts to partly compensate for the errors Introduced in such 
calculations by decoupling the extreme surface layer of the target 
from the recoil cascade, except for collisions which actually 
yield a sputtered particle. 

. This figure shows the fractionation predicted for a hypothetical 

target containing traces of mass 40 and 44 imbedded in a "background" 
matrix of atomic mass Mj- The dashed curve shows the fractionation 
calculated in the linear approximation (12), and the solid curve 
gives the exact fractionation calculated according to (2). The 
characteristic change over from negative to positive fractionation 
as a function of the background mass is clearly shown. Note that 
the non-linearity in the fractionation of the isotopes 40 and 44 is 
not given by the deviation f„rom a straight line of the curves shown 
here, since only the 6-values for a single pair are represented. 
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The curve shows the predicted fractionation In a pure Ca natal 
target containing terraatrlal abundances of the Isotopes. The 
deviation from linearity la due to higher order terns in the 
nass lncrenent c. The leading tern Is linear. 

The three curves In (a) and (b) summarize the fractionation 
pattern for three different choices of the relative scattering 
cross sections In an SiOg target. In (a), curve 1 corresponds 
to °o-</°Sl-0 “ °* B6 » curve 2 corresponds to ° 0 .(/o sl . 0 “ 1*16# 
and curve 3 corresponds to 0 0 -(/ 0 Si-0 " 1 ’°’ In curve 1 
corresponds to °si-Sl^°Si-0 “ curve 2 corresponds to 

°Si-Si /o Si-0 " °* 87 ’ * nd curve 5 corresponds to * sl _ sl /° st _ 0 " 

1 . 0 . 


In (a), (b), and (c) the lower panels show the fractionation 5 
for 0, Ca, and Ti, respectively, in perovsklte, CaTiO^. The 
upper panels in each case show the deviation from linearity as 
defined by (15). 

In (a), (b), and (c) the lower panels show the fractionation 6 
for 0, Ca, and Si, respectively, in anorthite, CaAlgSigOg. The 
upper panels in each case show the deviation from linearity as 
defined by (15). 

In (a), (b), (c), and (d) the lower panels show the fractionation 
5 for 0, Ca, Si, and Mg, respectively, for ackeroanite, Ca^MgSi^. 
The upper panels in each case show the deviation from linearity 
as defined by (IS). 


Fig. 9. In (a), (b), and (c) the lower panala show tha fractionation 6 
for 0, Si, and Mg, respectively, in enatatlta, MgSlOy Tha 
upper panala in aach caaa show the deviation fron linearity as 
defined by (15). 

Fig. 10. In (a) and (b) the lower panels show the fractionation & for S 
and Fe, respectively, in troilite, FeS. The upper panels in 
each case show the deviation fron linearity as defined by (15). 

Fig. 11. The Ca fractionation curves for CaFp (bottom) and Calg (top) 
are shown. The Insets depict the deviation of these curves 
from a straight line through the origin and through the 6(1*2,10) 
point. 
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